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Five and six dimensional static, spherically symmetric, asymptotically Euclidean 
black holes, are unstable under gravitational perturbations if their mass is lower 
^ \ than a critical value set by the string tension. The instability is due to the Gauss- 

Bonnet correction to Einstein's equations, and was found in a previous work on linear 
' stability of Einstein-Gauss-Bonnet black holes with constant curvature horizons in 

■ arbitrary dimensions. We study the unstable cases and calculate the values of the 

critical masses. The results are relevant to the issue of black hole production in high 
energy collisions. 
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I. INTRODUCTION 



The most conservative approach to gravity in higher dimensions is the one due to Lovelock 
[1], in which the LHS of Einstein's equation G ab + Ag ab = 8nG T ab is replaced with Q ab , the 
most general symmetric, divergency free rank (0, 2) tensor than can be constructed out of 
p the metric and its first two derivatives. Lovelock's tensor is 

j-h : icd-i)/2] 

Qab = Yl C " G(n ^ W 
n=0 



where D is the spacetime dimension, [z\ the highest integer satisfying [z] < z, and G^ n \b is 
obtained by making appropriate contractions on a tensor product of n copies of the Riemman 
tensor, contractions that trivially vanish if n > [(D — l)/2]. 

The first few GW^'s are the spacetime metric ab = g a b, Einstein's tensor G^\ b = 
R ab — \Rg a b-, and the Gauss-Bonnet tensor 

/^f(2) a p de p ca opcp da op CD o i p p a _a I p e/p cd a p d p c . p2\ 

Cx v 'fe — Kd, tide — ^tid ticb — ^tib tic + tXii-b ~ ~^9h {-tied tief ~ 4/£ c t^d TA j, 

(2) 

If D = 4, G^ n \b vanishes for all n > 1 and Lovelock theory reduces to Einstein theory with 
a cosmological constant Cq. Starting with D = 5, we may add the G^ 2 \ a term, and the 
resulting theory, usually referred to as Einstein-Gauss-Bonnet theory (EGB, for short), is 
the most general Lovelock theory in five and six dimensions: 

AG (0 ) b a + G(i) b a + ctG(2) b a = 8nG T b a , (3) 
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As is well known, EGB theory arises in the low energy limit of heterotic string theories [2, 3], 
a > being proportional to the inverse string tension, thus string related treatments of BHs 
in higher dimension should use the EGB equations. Spherically symmetric, asymptotically 
Euclidean vacuum black hole solutions of the EGB equations (3) with A = are well known 
since the eighties [3-5] . They are given by 



ds 2 = -f(r)dt 2 + f{r)- l dr 2 + r%c£rW, 
gijdx l dx j the line element of S n , n = D — 2, and 



(4) 



f(r) = 1 + 



4a/i (n - 1) in - 2) 



1 + 



nr 



n+l 



a(n - l)(n - 2) 

/i above is an integration constant, related to the mass M of the black hole through [6, 7] 



(5) 



M 



// 



8nG 



r(H±i) 



(6) 



A n being the area of the n— sphere. For positive \i and a, the case we are interested in, 
there is a single horizon located at the only positive root of (note the missing factor of 
1/4 in 



/i 



nr 



(n-3) 



[a(n- l)(n-2) + 2r 2 ] , 



then 



M 



nr h 



(n-3) A 



32ttG 



[a(n-l)(n-2) + 2r h 2 ] . 



The temperature and entropy of the black hole (4)- (5) are [6] 



T = 

S = 



n 



AG 



1 + 



2r h 2 + a(n - 2){n - 3) 
r h 2 + a(n- l)(n-2) 
an(n — 1) 
r h 2 



(7) 
(8) 

(9) 
(10) 



The specific heat can be obtained from (8) and (9) using 



Introducing f h :- 



C = - 



C 



we obtain 



DM 
W 



dM 



dr h J \dr h 



dT 



(11) 



nA n a n/2 
AG 



jn-2/ -2 



f 2 h + (n- \){n - 2)) 2 (2f 2 + (n - 2)(n - 3)) 
2r\ + (n-2)(n- 7)f 2 + (n - \){n - 2) 2 (n - 3) 



(12) 



Note that (4)- (5) reduces to the n + 2 dimensional Schwarzschild-Tangherlini [9] black hole 
of Einstein's theory in the a — > limit, since 



(13) 



3 



The thermodynamic functions (9)-(ll) reduce to the Schwarzschild-Tangherlini (ST) ones 
in the a — > limit. Note, however, that some solutions to the EGB equations are found 
to diverge as a — > 0, an example being the solution (4)-(5) with a plus sign in front of the 
square root in (5). Other crucial issues strongly depend on a being nonzero (we will restrict 
to a > 0, as in string theory). Consider first five dimensional (n = 3) EGB black holes. From 
(7) follows that there is a minimum mass /x = |a for black hole formation, otherwise, (4)-(5) 
has a naked singularity. This does not happen for five dimensional Schwarzschild-Tangherlini 
(ST) black holes. The temperature 

T = Th 
5D 27r(ri + 2a) 

goes to infinity as r h — > + (/i — > + ) for ST holes, whereas it tends to zero as r ft — > + 
(/i — > ) in the EGB case. The specific heat is always negative in the ST case, whereas it 
has a pole in the EGB case at r h = \[2a, with C > for r h < \[2a, and C < for r h > \[2a, 
i.e, small five dimensional EGB black holes can be in equilibrium with a heat bath, contrary 
to what happens for ST holes. Six dimensional EGB black holes behave more like ST black 
holes, their temperature decreasing monotonically from infinity in the interval < < oo, 
and their specific heat being always negative. However, both five and six dimensional low 
mass EGB black holes were found to be unstable under (linear) gravitational perturbations 
[8, 10, 11], whereas all D > 4 ST black holes are well known to be stable under linear 
gravitational perturbations [12]. 

The purpose of this work is to find the values for the critical mass below which five and 
six dimensional EGB black holes become unstable under linear gravitational perturbations. 
The perturbation treatment in [8, 10, 11] is based in the decomposition in tensor, vector 
and scalar modes given in [13], which is a higher dimensional generalization of the axial 
an polar modes found in the Regge- Wheeler treatment of Schwarzschild perturbations [14]. 
The metric perturbation in the tensor modes are made from symmetric, divergency free 
tensor fields on S n satisfying D k D k T^ = —k T 2 Tij, Dj the covariant derivative on S n . 
Similarly, vector (scalar) mode perturbations are made from vector (scalar) fields satisfying 
D]~D k Ti = —ky 2 Ti (DkD k T = —ks 2 T). A detailed exposition of the construction of these 
modes can be found in [13]. The spectrum of the Laplacian acting on divergency free, rank 
p symmetric tensors on S n is [15] 

k p 2 = £(£ + n -l)-p, £ = 0,1,2,... (14) 

ks, ky and fcr correspond to p — 0, 1 and 2 respectively. 



II. SCALAR MODE INSTABILITY OF FIVE DIMENSIONAL BLACK HOLES 

In five dimensions there is a single horizon located at (see (5)-(7)) 



r h = \/-fi-a, (15) 

as long as \x is greater than 3a/2, the minimum value required for black hole formation. It 
is convenient to adopt the dimensionless variables from Section 4a of reference [11] 

x := ra -1 / 2 , m := — , (16) 
a 
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then Xh = \J\ m ~ 1 an d the dimensionless tortoise coordinate 



x*{x) 



J2x h 



dx' 



f(x) = 1 + 



X 



1 - 



8m 
3x 



(17) 



extends from minus to plus infinity. 

Scalar perturbations in five dimensions (n = 3) of harmonic number ks 2 = £{£ + 2),£ = 
2,3,..., (the modes £ = 0, 1 are trivial [13]) are entirely described by a single function 
(f>(t, x) governed by an equation (7ik s + otd 2 /dt 2 ) = 0, which admits separation of variables 
<f>(t,x) = (f)(x)e u ' t , giving TCk s 4> — —ollo 2 ^) = aE<p, with (f> satisfying appropriate boundary 
conditions (reference [11], eqs. (61)-(66)). The "Hamiltonian" 



n 



d 2 



k s 



dx 



*2 



(18) 



can be constructed following Section V in [11]. A negative eigenvalue of Hk s -real u- implies 
that this mode grows exponentially with time, i.e., is unstable. Generic perturbations have 
projections on each harmonic (tensor, vector or scalar) mode. Since 5D BHs were found to 
be stable under tensor and vector perturbations [8, 10, 11], they will be unstable if and only 
if a ks is found such that the spectrum of TCk s is not positive. The boundary conditions 
defining the space of functions on which 7i ks acts determine its spectrum, L 2 (x*, dx*) being 
an appropriate function space for black hole spacetimes (see however, the discussion in [16] 
regarding nakedly singular spacetimes). The problem of stability is then entirely equivalent 
to the quantum mechanical problem of determining the sign of the lowest eigenvalue for 
each member of the family of Hamiltonians Hk s ,ks = \J £{£ + 2), £ = 2,3,... Our strategy 
to prove instability consists in showing that, if fi/a is small enough, then for sufficiently 
high ks, there exists a wave function with a negative expectation value of Hk s (numerical 
evidence of this fact was given in Section IVa of ref [11]). This implies that the ground state 
of TCk s nas negative energy, from where the instability follows. 

From the results in Section IV of [11], we find, after a long calculation, that, after introducing 



x := \/x h 2 + 1, y := \J x A + Ax 2 



the potential can be conveniently split as 



U, 



aVi 



k s 



ks 



k s 2 Qoc + qo + 



ks 2 q\ + <?2 



where V is a quartic polynomial in ks'- 

V = 2x 2 y 4 [(k s 2 - 3)y + 6x 2 )] 2 , 



(19) 



(20) 



(21) 
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and the q's do not depend on ks'- 
{x 4 - Ax 2 ) 



goo 



Qi 

q-2 



X 2y2 

(x 8 + 120x A x o 2 - 240x o 4 ) 
8x 2 y 3 

V° - 6x 8 + 200x o 2 a; 6 + 528x a 2 x 4 - 560x o V - 480x o 4 

8x 2 y 4 

24x 2 x D 2 [(24x 2 + 20x o 2 x 2 - 6a; 4 - x Q ) y - A8x Q 4 - 8x D 2 x 4 + x 8 ] 
72x 2 Q x 4 [(x 6 - (2x 2 - 6)x 4 - 20x 2 o x 2 - 2Ax 4 Q - 2Ax 2 ) y 
-x 8 + 2x 2 Q x 6 + Ax 2 Q x 4 + Attx 4 x 2 + 96x 4 ] 



(22) 



(23) 
(24) 

(25) 



Note that q^ is negative in the range \x\ < x c := \/2 (1 + Xh 2 ) 1 ^ 4 , and that < Xh < x c 
if and only if 3/2 < m < 9/2 + 3y/2. Suppose this is the case and let tp(x) be a real C°° 
function vanishing outside (xh,x c ), normalized such that 



ipip dx* = 



1 



dx 



Using tp as a test function, the expectation value of the kinetic piece of (18) is 



< -d 2 /dx* 2 >= - 
and that of the scalar potential is 



i/j f ^ dx* 



dx 



*2 



f 







,dijj 



— f— 

dx \ dx 



dx 



(26) 



< aVk s >= a 



f 

J — c 



i/jV ks i/jdx* 



rx c 

/ i) 2 U ks dx = k s 2 Qoo + Qo + Q(k s ) 



where 



and 



Qoo= ^ 2 qoodx < 
Qo= ip 2 q dx 



do not depend on ks, and 



Q(ks) = / 



,: 2 



ks 2 q\ + Q2 
V 



dx 



(27) 

(28) 
(29) 

(30) 



Note that the integrand in (30) converges uniformly to zero in the interval x G [a^xj as 
ks — > oo. This follows from the fact that T> is strictly positive in [x/j,a; c ] (see (21)) and is 
a quartic polynomial in ks- As a consequence, \vai ks _ 00 Q{ks) = and thus < Ti. ks > is 
negative for the given test function and large values of ks- Since the above construction is 
possible if 

3/2 < m < m c{5D) = 9/2 + 3^2 ~ 8.743, (31) 
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we conclude that, in this mass range, all (static, spherically symmetric, asymptotically 
Euclidean) 5D black holes have a high harmonic scalar instability. 

Although solving the quantum mechanical problem (18) analytically is out of consideration, 
in some cases we were able to spot the fundamental energy using a shooting algorithm to 
numerically integrate (18). This was done in the standard coordinate x (instead of x*), for 
which (18) reduces to an equation of the form <p" + P<f>' + Q(f) = with a regular singular 
point at the horizon. The first few terms of the Frobenius series around the horizon were 
used to generate appropriate initial conditions for the shooting algorithm. As an example, 
we exhibit in fig.l the scalar potential vs. x, together with the ground state wave function 
corresponding to m — 1.7, t — 2. We also remark that no bound state was found for 
m > m c{5D) . 
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FIG. 1: The scalar potential (arbitrary scale) in five dimensions for £ = 2 and fi/a = 1.69 is shown 
together with the ground state wave function and energy, found numerically using a shooting 
algorithm. The origin of x* was chosen as in (17), x* = for x = Ix^. 



III. TENSOR MODE INSTABILITY OF SIX DIMENSIONAL BLACK HOLES 

As in [8], we find it convenient to introduce dimensionless variables 

m := fia~ 3/2 , x := r/{fia) 1/5 = ra^m- 1 / 5 , dx*/dx := 1/f (32) 

The spectrum of the Laplacian on symmetric divergency free tensors on S* 4 is kr 2 = l{l + 
3) — 2, t G Z, only t > 1 tensors being required to construct non trivial tensor perturbations 
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of 6D black holes. These perturbations are entirely described by a single function 4>(t,x) 
governed by an equation that, after separation of variables <fi(t,x) = <f)(x)e UJt , assumes the 
form Hk T ?p = —am 2 ^uj 2 <\) = am 2 ^ J E(j) ([8], equation (16)), with "Hamiltonian" 



a 



= ~ft^2 + V **> ( 33 ) 

Vk T being the RHS of eq. (18) in [8]. From [8] we can readily construct the potential, the 
result is 

U := = (A; T 2 + 2)m-i U a + m"* U x + U 2 - U 3 (34) 

where the U/s depend only on x: 

2(g 5 + 6 ) 2 - 75 

2x 2 (x 5 + l)(x 5 + 6) 1 ' 
8 x 20 + 72 x 15 + 1218 x 10 + 1752 x 5 - 27 

1 ~ 4x 2 (x 5 + l) 2 (x 5 + 6) 2 ( ' 

24 x 20 + 336 x 15 + 2414 x 10 + 2916 x 5 + 189 

2 ~ 24(x 5 + l) 2 (x 5 + 6) 2 (3 ' 

24 x 20 + 216 x 15 + 1154 x 10 + 1506x 5 - 81 

Us — , = (38) 

24x 5 (x 5 + l) 2 (x 5 + 6) y/l + 6/x 5 

Let x c = (a/75/2 — 6) 1 / 5 ~ 0.658 be the only positive root of U , note that U < for 
< x < x c . The x coordinate of the horizon is 



z 

Xh 



2 _ A , x 1/3 



(2m + 2 v / 16 + m 2 ) . (39) 



2ms z 

Xh is a monotone increasing function of m, and Xh = x c at m = m c (6D) given by 

72V6(5V / 6-12) 3/2 
(l2 - 5A/6 + 6 x /4 a/5 A/5A/6 - 12) 



m c(6D) = ^ ; ^ 7.965 (40) 



If m < m c ( 6 £)), then x^ < x c and we can take a test function supported in (x/j,x c ), so that 
the expectation value of the Uo piece of the potential is negative. Note from (34) that this 
term is proportional to the harmonic 1(1 + 3), and no other term of 7ik T depends on £, thus 
the expectation value of Hk T for such a test function will be negative for sufficiently high 
harmonic number. We conclude that 6D BHs are unstable if m < m c above. 
Now we prove stability for m > m c ^ D y. U Q ,U\ (and U2) are positive if x > x c , whereas 
U2 — U3 > if x > x' c ~ 1.176. Since Xh = x' c for m = m! ~ 48.927, stability will follow if 
we prove that U > for £ = 2,3..., m c ( 6D ) < m < m' and x > Xh given in (39). A lower 
bound for U in this region of parameter space is given by the minimum of the single variable 
function U L := (10U + £/i)50~ 2/5 + U 2 - U 3 in the interval x e (x c , 00). After some work 
Ul can be seen to be positive in this interval, thus proving stability. We conclude that 6D 
BHs are linearly unstable if and only if fi/a 3 ^ 2 = m < m c (6D)- 

Fig. 2 exhibits the potential and fundamental state (found numerically) corresponding to 
i = 2, fia~ 3 / 2 = 1.85. 
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FIG. 2: The tensor potential (arbitrary scale) in six dimensions for £ = 2 and /j,a~ 3 / 2 = 1.85 is 
shown together with the ground state wave function and energy, found numerically using a shooting 
algorithm. The origin x* = corresponds to x = 2xh- 

IV. CONCLUSIONS 

Gauss-Bonnet corrections to Einstein's equations in higher dimensions have been consid- 
ered in many different models, and naturally arise in the low energy effective action of certain 
string theories. However, their effects on black hole formation has long been disregarded. 
The instability found in [8, 10, 11] and this paper implies that the simplest EGB black holes 
(asymptotically Euclidean, static, spherically symmetric), which are the closest analogue of 
Schwarzschild black holes, cannot actually be formed in five space time dimensions if their 
mass parameter fi (see (5)-(6)) is less than ~ 8.743a. The Gauss-Bonnet term also prevents 
the formation of these black holes in six dimensions unless fx is greater than ~ 7.965a 3//2 . The 
implications of these figures depend on the context where (4)-(5) is used. As an example, 
the n-dimensional EGB black hole (4)- (5) is an approximate EGB solution if we periodically 
identify one of the asymptotically Euclidean coordinates with a period much larger than the 
horizon radius, and our perturbative analysis should be valid in this setting. The large extra 
dimensions scenario (suitable only for D > 6, [17]) is of interest because it allows a to be in 
the TeV scale [17, 18], and so mini black holes could be produced in high energy collisions 
and be eventually detected at LHC. In view of our results, the probability of these events 
may be severely limited due to low mass black hole instabilities. As far as we know, this 
fact has not been taken into account in previous calculations on black hole production rates 
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in high energy collisions. In theories where the EGB equations simply arise as a low energy 
effective theory of some quantum gravity model, a is of the order of the Planck scale and 
the bounds we obtained for small black hole masses are much more stringent. 
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